Critical Behavior of the Two-Dimensional Random Quantum Ising Ferromagnet 
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We study the quantum phase transition in the two-dimensional random Ising model in a transverse 
field by Monte Carlo simulations. We find results similar to those known analytically in one- 
dimension: the dynamical exponent is infinite and, at the critical point, the typical correlation 
function decays with a stretched exponential dependence on distance. Away from the critical point, 
there may be different exponents for the divergence of the average and typical correlation lengths, 
again as in one-dimension, but the evidence for this is less strong. 



I. INTRODUCTION 

Though classical phase transitions occurring at finite 
temperature are very well understood, our knowledge of 
quantum transitions at T = is relatively poor, at least 
for systems with quenched disorder. There is, however, 
considerable interest in these systems since they (i) ex- 
hibit ncwjiniversality classes, and (ii) display "Griffiths- 
McCoy"ErB singularities even away from the critical point, 
due to rare regions with stronger than average interac- 
tions. 

Just as the simplest model with a classical phase tran- 
sition is the regular Ising model, the simplest random 
model with a quantum transition is arguably the Ising 
model in a transverse field whose Hamiltonian is given 
by 



{id) 



fX,- cr. 



(1) 



Here the {cf } are Pauli spin matrices, and the nearest 
neighbor interactions and transverse fields hi are both 
independent random variables. Naturally this model has 
been quite extenaivelir studied and many analytical re- 
sults are availableO~Ero for the case of the |One-dimensional 
chain. Among the surprising predictiona3 for the critical 
behavior in 1-d are: 

1. The dynamic critical exponent, z, is infinite. In- 
stead of a characteristic time scale varying as a 
power of a characteristic length scale ^ according 
to ~ one has instead an exponential relation 



£^r ^ exp (const. ^^), 



(2) 



where n 
scaling. 



1/2. This is called activated dynamical 



Distributions of the equal-time a^-a^ correlations, 
are very broad. As a result average and typicalu 
correlations behave rather differently, since the av- 
erage is dominated by a few rare (and hence atyp- 
ical) points. At the critical point, for example, 
the average correlation function falls off with a 



X = 2 — (p ~ 0.38 with the golden mean, whereas 
the typical value falls off much faster, as a stretched 
exponential 



with 



Ctypir) ^ exp(-const. r"), 
1/2. 



(3) 



Away from the critical point the average and typi- 
cal correlation lengths both diverge but with differ- 
ent exponents, i.e. ^av 5'^""] ^typ 5'^*-^^, where 
5 is the deviation from criticality and 



= 2: 



Vtyp = 1- 



(4) 



power of the distance r as Ca.v{r) 



where 



Unfortunately, the analytical approach is only valid 
in 1-d, and very little is known about the critical be- 
havior in higher dimensional. An important question is 
whether these striking analytical results are a special fea- 
ture of 1-d, or whether they are valid more generally. 
Here, we make a first step in investigating this ques- 
tion by performing large-scale Monte Carlo simulations 
of the random transverse field Ising ferromagnet in two 
dimensions. Because the ferromagnet has no frustratig 
we are able to use highly efficient cluster algorithmsQ' 
which considerably reduce critical slowing down. Conse- 
quently, we are able to study larger s izea and get much 
better statistics than in earlier workE3llil on quantum 
spin glasses. Furthctmore, cluster algorithms have an 
improved estimatoill3 for correlation functions, which 
greatly reduces the statistical errors when the mean is 
small, so we are able to probe the distribution of corre- 
lation functions out to relatively large distances. 

Our main conclusion is that the behavior of the 2-d 
system at the critical point is very similar to that of 
the 1-d case, see points |l] and ^ above. There are some 
difficulties in interpreting our data away from the critical 
point but, nonetheless, here too there is some evidence 
for behavior similar to that in 1-d, see point ^ above. 
Further details and additional results will be presented 
elsewherellj. _ 

In related work, Rieger and KawashimallJ have studied 
the same model but with continuous imaginary time, con- 
centrating on Griffiths-McCoy singularities in the para- 
magnetic phase, which are characterized by a continu- 
ously varying dynamical exponent, z{5). Kawashima and 
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Rieger find that z{6) tends to infinity as the critical point 
is approached {6r^ 0), again as in 1-cH. 

As is standardta, we represent the d-dimensional quan- 
tum Hamiltonian in Eq. by an effective classical 
"Hamiltonian " in (c?-|-l)-dimensions, where the extra di- 
mension, imaginary time, is of size j3 = 1/T and is di- 
vided up into Lr intervals each of width At — (3/Lr. 
The effective classical Hamiltonian is given by 

PciHci^^ if,,5,(T)5,(T) -^i^,5,(r)^,(r') (5) 

where, r' = r + At, Si{T) = ±1, Kij = ArJij, 
exp(— 2if,;) = tanh(AT/ii), and Pel = ^/Td where Td 
is an effective "classical" temperature (not equal to the 
real temperature which is the inverse of the size in the 
time direction). Note that the disorder is quenched and 
so the interactions are independent of t. Correlations of 
the classical Ising spins, S'i(T), correspond to a^-a^ cor- 
relations in the original quantum model. These are the 
correlations that we are interested in here. 

In order to capture the random quantum critical be- 
havior in the intermediate size systems that we can sim- 
ulate, we wish the disorder to be effectively quite strong. 
In particular, we would like some of the fields to be much 
stronger than the bonds in their vicinity and vice-versa. 
This is captured by having distributions for both the 
fields and interactions with a finite weight at the origin. 
This is satisfied by a constant distribution for the Kij and 
a distribution for the Ki which falls off exponentially at 
large values. Furthermore, the critical behavior does not 
depend on short distance physics and is hence indepen- 
dent of At, so for convenience we will take At = 1. A 
Hamiltonian which encapsulates these features in a way 
convenient for simulations is 

Hcl = - J2 J^3S^{r)S,{T) - ^ J,5.(t)5,(t + 1), (6) 

where t runs over integer values, 1 < t < L,-, and the 
distributions of the interactions are given by 



7r(J) 
PiJ) 



1 for < J < 1, 
otherwise, 

2exp(-2j) for J^>0, 
forJ<0. 



(7) 



Disorder averages are denoted by [• • -Jav and Monte Carlo 
averages for a single sample are denoted by (•■•). One 
controls the strength of the fluctuations, and hence tunes 
through the transition, by varying the classical tempera- 
ture. Tel. 

To assess the accuracy of the method, we not only 
performed simulations in 2-d, but also carried out sim- 
ulations for the 1-d model for a similar range of sizes 
and compared the results with the analytic predictions 
discussed above. The size of the system in the space di- 
rection denoted by L, so there are N = L'^L^- sites in the 



sample. We studied sizes up to L = 48 and Lr = 2048 us- 
ing the WolflEI cluster algorithm. We found that no more 
than 100 sweeps were required for equilibration, even for 
the largest lattices. At least 1000 samples were averaged 
over. 

The first task is to locate the critical point, for 
which we use ajjnethod suggested by D. Huse (private 
communication)L3lid. One computes the Binder ratio 



<?av - 2 



(M2p 



(8) 



where M = X)i t (''')' "^tiich (assuming, for now, that 
z is finite) has the finite-size scaling form 
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(9) 



For fixed L, gav ts^s a peak as a function of Lr- At the 
critical point, TcEj, the peak height is independent of L 
and the values of Lr at the maximum, L^'^^, vary as L^. 
Furthermore, a plot of gav against Lr/L™^'^ at the critical 
point, which has the advantage of not needing a value 
for z, should collapse the data. We see in Figs. I and I 
that this does not happen either in 1-d or 2-d. Rather 
the curves clearly become broader for larger sizes. This 
is easy to understand since we know that for 1-d z is 
infinite and it is the log of the characteristic time which 
scales with a power of the length scale. This suggests 
that the scaling variable should be h\Lr / hiL™^^ with 
InL™^'^ ~ L^, say. The inset to Fig. ^ shows that this 
works moderately well for 1-d (though not perfectly for 
this range of sizes) with the expected value yit = 1/2. In 
2-d, the data collapse for sizes L > 12, shown in the 
inset to Fig. |[ is quite good for = 0.42, and not quite 
so good with the 1-d value, /i — 1/2, though we would 
not claim that /i = 1/2 is ruled out. The close similarity 
of the data for 1-d and 2-d, particularly the broadening 
of the data in the main part of the figures, suggests that 
z is infinite also in 2-d. 

Next we consider the behavior of the equal time 
correlations at the critical point using the improved 
estimatorll^ _ Figs. ^ and ^ show data for the average 
and typicalQ correlations for spins separated by r = L/2 
in 1-d and r — {L/2, 0) in 2-d. For each value of L, we 
took Lr such that g^cv is close to the pealJlJ, see Figs. |l| 
and 1^. According to finite size scaling, the dependence 
on L for a finite system should be the same as the de- 
pendence on r in a bulk system. Fig. |^ shows that in 
1-d, the average correlation falls off with a power law 
while the typical valuCpdecays with an exp(— const. L^/^) 
variation, as expectedcl. Note, though, that the expo- 
nent for the average, —A, obtained by fitting the data- 
gives A = 0.50, somewhat different from the expectedcl 
value of 0.38. We will discuss this discrepancy further 
below. The data in Fig. ^ shows that the behavior in 2-d 
is quite similar: the average falls off with a power law, 
with A about 1.5, while the typical value falls off faster 
than a power law, (because of the downward curvature) 
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FIG. 1. A plot of the data for the average Binder ra- 
tio, 5av, at the critical temperature Td = Tc — 0.98 for the 
one- dimensional model. In the main figure the horizontal axis 
is L^/L™"" where L™"^ is the value of Lr at the peak. Note 
that the curves do not scale but rather get broader for larger 
sizes, indicating activated scaling, 2; = oo, (which is known 
from the analytical workB). The inset shows a plot appro- 
priate for activated scaling in which the horizontal axis is 
In L-r/ In L^"', with In L?'^" ~ L^/^ , which is also known from 
the analytical work. Scaling works quite well though there are 
still some small deviations for this range of sizes, even in the 
large Lt limit. 



consistent with a stretched exponential behavior of the 
form exp(— const. L'"), with a ~ 1/3. The statistical er- 
rors (shown) are generally smaller than the size of the 
points, so the downward curvature is statistically signif- 
icant. The value of a seems to be different from the 1-d 
result of 1/2 (the data is curved on the appropriate plot) 
but we are not able to determine it with great precision. 

Finally we briefly describe our data away from the crit- 
ical point, keeping keep the same values of L and L^. as 
at the critical point. In 1-d, (7av scales well with the ex- 
pected value I'av — 2 and the corresponding plot in 2-d 
also works well, this time with — 1.5. Plots of the 
typical correlation function (divided by the value at crit- 
icality) give an exponent close to the expected value of 
'^typ = 1 in 1-d and around 0.7 in 2-d, i.e. in each case 
much smaller than the exponent for the average correla- 
tion length obtained from pav However, attempts to de- 
termine I'av directly by scaling data for the average corre- 
lation function away from criticality were less successful. 
In 1-d, the fit, which works quitCr^vell, gives I'av — 1.2, 
much less than the expected valueH of 2, while in 2-d we 
found Vav — 0.75, not significantly different from vtyp- 

This discrepancy for results from average correlation 
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St/ t 
but for two dimensions, with 

Tci — Tc — 2.45. The conclusion is similar, i.e. the broadening 

of the curves for large L in the main plot indicates activated 

dynamical scaling, i.e. z = 00. The data for L = 48 is slightly 

high which may indicate that the true value of Tc is a little 

higher. In the inset, the data for L > 12 is seen to scale quite 

well with the same form InLr/I/'' known to be exact in 1-d, 

but the value of fi = 0.42 (shown) works better than the 1-d 

value of = 1/2. 



function may be due to difficulties in taking enough sam- 
ples to probe the rare regions which dominate the aver- 
age. This may also be why the exponent for the decay of 
the average correlation at criticality in 1-d was not found 
accurately. Another concern with the data away from the 
critical point is that we do not understand the way that 
the critical singularities go over to Griffiths-Mdjity sin- 
gularities, which have a continuously variableclolla z{S). 
Perhaps the latter give big corrections to scaling, thus 
possibly making the size of the critical region quite small. 

To conclude, we have found a strong similarity between 
the critical behavior of the random transverse field Ising 
model in one and two dimensions. The evidence that 
z ~ 00 and that the typical correlations at the criti- 
cal point decay with a stretched exponential function of 
distance are quite strong, but the evidence for different 
exponents for the average and typical correlations in 2-d 
is weaker because of discrepancies involving the data for 
the average correlation function. 

After this work was completed, we heard that S.-C. 
Mau, O. Motrunich and D. A. Huse (private commu- 
nication) have implemented numerically for d > 1 the 
renormalization group approach used in Ref. ||, finding a 
flow to the infinite disorder critical fixed point, just as in 
d=l. 
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FIG. 3. The main figure shows the average and typical^ 
correlations between spins L/2, apart at the critical point in 
one dimension. As expectedtl, the average decays with a 
power law (though the fit to the data gives a slope, —A, with 
A = 0.50 compared with the analytical value of 0.38), while 
the typical value decays faster than a power, as shown by the 
downward curvature in the plot. The inset shows that the 
typical correlation function decays as the stretched exponen- 
tial, Ctyp(L/2) ~ exp(— const. L^^^), as expectedQ. 
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